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Abstract 

We derive a theoretically allowed domain for the charge radius p and curvature 
c of the Isgur-Wise function describing the decay A& — ► AJui- Our method uses 
crossing symmetry, dispersion relations and analyticity in the context of the Heavy 
Quark Effective Theory but is independent of the specifics of any given model. 
The experimentally determined values of the T masses have been used as input 
information. The results are of interest for testing different models employed to 
calculate the heavy baryon formfactors which are used for the extraction of \V c b\ 
from the experimental data. 
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1. INTRODUCTION 

During the last few years, considerable progress has been made in the development of 

methodologies for the extraction of \V c b\ [1-6] from both exclusive and inclusive semilep- 
tonic B meson decays, where the leptonic current is clearly separated from the matrix 
elements of the hadronic current. The CKM matrix element V c b has important implica- 
tions for the investigation of rare decays and CP violation. It is therefore imperative to 
have as many independent and accurate determinations of \V c b\ as possible. There are 
however some theoretical limitations in the determination of \V c b\ from inclusive semilep- 
tonic decays due to issues related to validity of quark-hadron duality near the kinematic 
endpoint region 0. On the other hand, the exclusive decays must be described in terms 
of a number of undetermined nonperturbative formfactors that contain the physics of the 
hadronisation process. 

From the theoretical point of view, in the limit where the heavy quark mass tends 
to infinity, the analysis of heavy baryonic semileptonic decays are comparable to that of 
semileptonic decays of heavy mesons. This is due to the fact that heavy quark symmetry 
in the above limit predicts a single universal formfactor for both these types of decays 
known as the Isgur-Wise (IW) function ||. It is therefore interesting to estimate V c b from 
the decays of heavy baryons. However, experimental data on the decay of heavy baryons 
is still sparse compared to the corresponding meson decays. Semileptonic decays of the 
Ab baryon have up to now only been observed at LEP || [Kj. It is however expected 
that in the near future more data would be available from the LEP as well as from the 
the forthcoming B factory. Once these results are available, the theoretical value for the 
semileptonic decay width of A& can be used to gain an independent determination of \V c b\. 

In semileptonic decays of both heavy mesons and baryons, the heavy quark symmetry 
predicts the value of the corresponding IW functions at zero recoil || fTI] , though the 
shapes of the IW functions are left unspecified and require detailed knowledge of the 
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non-perturbative strong interaction physics. Unfortunately, the differential decay width 
vanishes at zero recoil and hence the known normalization of this function at that 
point is not enough to extract V c b from the data. To overcome this difficulty various 
models are used to determine the functional form of the IW function and is extracted 
therefrom. However, the observed variations of \V c b\ obtained from different models turn 
out to be larger than the statistical and systematic uncertainties in the experiments. 
From the experimentally measured lepton invariant mass spectrum one may determine 
in a model independent way by extrapolating to the kinematical endpoint of maximal 
momentum transfer to the leptons corresponding to the zero recoil point. However, this 
extrapolation requires a large amount of data, very close to q^ax-i which is difficult to assess 
experimentally. Thus measurement of \V c b\ requires a parametrization of the IW function. 
It is therefore necessary to put constraints as far as possible in a model-independent 
manner on this function. 

In this paper we present an analysis of the model independent constraints on heavy 
baryonic I W function arising in the description of semileptonic A& — > A c lv\ decay. This 
is accomplished by first obtaining constraints on the charge radius and the convexity of 
the I W function for the elastic A& formfactors using the techniques of dispersion relations 
and analyticity. Next, heavy quark symmetry is used to relate this IW function to the 
formfactors appearing in the semileptonic A& — > A c Wi decay. The dispersion theoretic 
method of extracting information on amplitudes is quite old |12| , [13], [14| and has been 
applied to the study of semileptonic decays of light mesons in a more contemporary 
language |L5], Its application to heavy quark systems has received much attention in 
more recent years [7, 17 - 24]. In Section 2 we give a brief account of the general formalism 
for the decay of heavy baryons. The method for the estimation of model independent 
bounds on the baryonic IW function is described in Section 3. The numerical results 
are presented in Section 4, and Section 5 concludes the paper with some discussions and 
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future outlook. 



2. GENERAL FORMALISM OF HEAVY BARYON DECAY 

In this section we give a brief account of the semileptonic decay of A& baryon using 

heavy quark symmetry ||26|| . 

The hadronic matrix element for the A b — > A c Wi decay is parametrized by the following 
general decomposition consistent with Lorentz invariance. 

(A c (p')\J,\A b (p)) = ^{[Mq^ + tUq 2 )^^ + f 3 (q 2 )q,} 

-[#i0 2 )7m + m{q 2 )(r^q u + 93(q 2 )q^]l5}uA b} (2.1) 

where J M = (V^ — A^) = 07^(1 — 7s)& and q — (p — p') is the momentum transferred to 
the leptons. Here u Ab and u Ac are the Dirac spinors of and A c respectively. The f's 
and g's are the formfactors for the independent Lorentz covariants corresponding to the 
hadronic matrix elements of the vector and axial vector currents respectively. 

When both baryons are heavy, it is convenient to parametrize the matrix element in 



terms of the four velocities v and v' [26[] : 

(A c (0|J M |A 6 (?;)) = u Ac {F l {u) lfl + F 2 {uj)v fl + F 3 {uj)v , II 

-[Gi(o;)7 (U + G 2 (lo)v^ + G 3 (u)v' ^ 5 }u Ab , (2.2) 

where the formfactors [written with corresponding uppercase F's and G's] are functions 
of u) — v ■ v', the only non-trivial Lorentz scalar formed from v and v'. 

In terms of these formfactors, the differential decay width in the zero mass approxi- 
mation for charged lepton is given by 
dT IVc^Gf 2 ^ 



dq 2 96n 3 m Ab 

where 



-{t-^mi 2 + \H V \ 2 } + t + [2t|G!| 2 + \H A \ 2 ]}, (2.3) 



Hy(q 2 ) = (m Ab + m^F, + ^( A. + -A), (2.4) 

2 m Ab m Ac 



H A (q 2 ) = (m Aft - mjjG, - U^- + -^_), (2.5) 

with t = q 2 , t± = (m\ b ± m Ac ) 2 — t, and k = \\J t - ± ^- 

In the infinite quark mass limit, the spin-flavour symmetry for the b and c quark 
relates all the formfactors to a single universal function which is the IW function. 

F 1 {u)=G 1 (u) = £(u), 

and 

F 2 {uj) = F 3 (uj) = G 2 {uj) = G 3 (u) = 0. 

The IW function so defined is independent of heavy quark masses and normalized at zero 
recoil (v = v') i.e. ^(1) = I. 

The behaviour of the IW function close to zero recoil is of particular interest and it is 
customary to parametrize this behaviour through a Taylor's series expansion about the 
zero recoil point namely, 

= 1 - p 2 (cu - 1) + c(u - l) 2 + 0[(u - l) 3 ], (2.6) 

where p and c are referred to as the charge radius and the convexity parameter respectively. 
The kinematic region, accessible to the semileptonic decay lies in such a small domain 
(w — 1 to 1.43) that a precise knowledge of the charge radius could basically determine 
the IW function in the physical region and the convexity parameter gives the corrections 
to the simple linear dependence. In the following section we derive bounds on p and c 
and thereby provide model independent constraints on the IW function. 

3. DERIVATION OF THE BOUNDS 

We consider the elastic formfactors of the baryon through the matrix elements of 

the flavour- conserving current = bq^b 

(A b (p')\V,\A b (p)) = « Ab {FfM 7M + F 2 e (uj)(v + v'),}u Ab . (3.1) 
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In the heavy quark limit, the b-number conserving elastic formfactors are related to the 
same IW function which enters in the expression for the differential decay rate 

A b — > A c lu. The short distance and the finite mass corrections are in this case much 
smaller than that for the formfactors involving flavour changing currents, corresponding 
to the decay A5 — > A C W. 

The method is based on the properties of the two-point function H(q 2 ) defined as 

W u (q 2 ) = i J rf 4 xe ig - x (0|T(^(x)\/ i/ (0)|0) 

= (gV-gV^)n(g 2 ), (3.2) 

with = Jry^b. The conservation of the current V 1 leads to the transverse nature of 
the two-point function. In QCD, the asymptotic nature of the two-point function IT(g 2 ) 
is such that it satisfies a once subtracted dispersion relation 

2 dn(g 2 ) , 1 roo ImII(t) 

The absorptive part ImII(t) is obtained from the unitarity relation 

W - g 2 Olmn(g 2 ) = i^^ r (27r) 4 ( 5( 4 )(g-j 5r )(0|(^(0)|r)(r|^(0)|0), (3.4) 

z r 

where the summation is extended over all possible intermediate hadronic states with 
the quantum numbers of the current and is weighted by the corresponding phase 
space factor dfi r . A judicious choice of /x and v (/i = v) makes this a sum of positive 
definite terms and thus one can obtain a strict inequality by retaining only the term 
with intermediate T = A^Afe state. From crossing symmetry, the A&A& — > vacuum matrix 
element is described by the same set of elastic formfactors, but relevant to the pair- 
production region (4m^ < q 2 < oo) instead of the elastic region (q 2 < 0). This gives rise 
to an integral inequality 
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where 

H*(q 2 ) = 2m Ab F 1 E + (Am 2 Ab - t)£-. (3.6) 

In the case of mesons, inclusion of other two-particle intermediate states related by spin- 
flavour symmetry were made to saturate the dispersion relation in order to find the rele- 
vant bound. For A& baryon decay, there are no spin symmetric partners that contribute to 
the dispersion relation bound. However, the presence of more than one helicity amplitude 
fulfills in effect the same role. 

For Q 2 far from the resonance region (2rri} ) AQCD 4m^ + Q 2 ), the two-point function 
can be computed reliably from perturbative QCD. For large b quark mass, it is sufficient 
to take Q 2 = 0. The one loop expression for x(0) is 



3 r 1 x 2 (l-x 



2^ 



Perturbative a s corrections to this result are negligibly small at the physical m\ b scale. 
The non-perturbative corrections |27 are included by expressing the two-point function as 
an operator product expansion and incorporating the leading nonperturbative gluonic and 
quark-antiquark vacuum condensate terms (G 2 ) and {qq). But these terms are suppressed 
by the fourth powers of a large mass scale for dimensional reasons. 

The analysis simplifies by using a conformal transformation to map the full complex 
t plane onto the unit disc in the complex z plane, 

l + z 



\ 



1--^-. (3.8) 



The transformation maps the cut [Am\ b , oo] on the unit z circle and the rest of the t plane 
into the open unit z disc. In terms of the new variables z the inequality eq. (3.5) gets 
translated to 

-L jf ^.[\^{z)Ff (z)\ 2 + |^(z)fl?(z)| 2 ] < 1, (3.9) 
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where the functions <f> F E and <p H E are defined as 



and 



= + z){l-zf'\ (3.10) 



= ^^(l + z)(l-zf 2 m A \ (3.11) 



These functions are analytic and nonzero inside the unit disc \z\ < 1. Their moduli 
squared on the boundary are equal to the positive weights appearing in the integrals, 
multiplied by the Jacobian | % | of the conformal transformation. To translate the bound 
into one that is valid in the semileptonic region, we need a function which is analytic inside 
the unit disc. The functions 0, have no poles, branch cuts or zeros in the interior of the 
unit circle \z\ < 1, but the formfactors F E (q 2 ) and Hy(q 2 ) have poles due to the existence 
of spin-one bb states (J p = 1~) below the A^A^ threshold. They also have branch cuts 
originating from non-resonant continuum contributions with the invariant masses below 



AfeAfe threshold, though their effects may be neglected |2T], [23j . The relevant states giving 
rise to poles in this context are the 3 Si and 3 Di bottomonium states, i.e. T(nS) and 
T(nD), respectively. From the experimental values of the masses of the T resonances, 
the position of the poles are obtained at the points: 

z x = -0.29, z 2 = -0.37, z 3 = -0.43, z 4 = -0.48, z 5 = -0.58, z 6 = -0.65. (3.12) 

The residues of the formfactors at these poles are however unknown. To make up for this 
lack of information, one introduces a product of the functions of the form (z — Zi)/(l — Ziz), 



known as Blaschke factors [28 



This function P(z) has the virtue that it is analytic on the unit disc \zi\ < 1, and elim- 
inates, through multiplication, the poles of F® at each z = zj. Since each term in P(z) 



is contrived to be unimodular on the boundary, the inequality remains unchanged. Then 
the inequality (3.9) becomes 

I -U FiE (z)P(z)F 1 E (z)\ 2 + |0^)P(z)P^(,)| 2 ] < 1. (3.14) 

It is important to note that up to now the formfactors are treated as completely unknown 
functions, because heavy quark symmetry is quite inapplicable in the timelike domain 
near the pair-production threshold , though is only valid in the neighbourhood of the 
zero recoil point. 

Both the functions <j>i{z) and P(z)F E (z) are now analytic in the unit disc. We can 
hence apply the well known results of interpolation theory |28j for vector-valued analytic 
functions to obtain the bounds on the formfactors at points inside the unit circle. Accord- 
ingly, we apply the inequality of the Schur-Caratheodory type [28] at the origin retaining 
terms up to the first derivative with respect to z 

0^(O)P 2 (O)F^(O) + 0^(O)P 2 (O)^(O) + (0 FiE PFf)l 2 (O) 

+ (<j> H zPH E y?(0) < 1, (3.15) 

where the prime denotes differentiation with respect to z. Through the conformal map- 
ping, the point to = 1 is mapped to z = 0. We can hence use the prediction of heavy 
quark symmetry for the formfactors at the point of zero recoil, which leads to 

Ff(z = 0) = £{z = 0) = 1, H E (z = 0) = 2m Kb i(z = 0) = 2m Ab . (3.16) 
The charge radius and convexity of the IW function are defined by 

? 2 = -[^U, -=#L„ (3.17) 



the derivatives here being with respect to uj. They are related to derivatives with respect 
to z through 

[^U = -8, 2 , [^U = I28c-3V. (3.18) 
8 



From Eqs. (3.15 - 3.18), the upper and lower bounds on the charge radius are readily 
obtained to read 

b-Vb 2 - 4ad 2 b + Vb 2 - Aad . . 
< p < (3.19) 



with 



o = 64P 2 (0)[^ + ^], (3.20) 
6 = 8P(O)[0 Ff (O)(0 Ff P)'M + 0^(O)(0^P);(O)], (3.21) 



and 



d = P\^W f e + r H *} + (c/> Ff Py z \0) + (0h^); 2 (O) - 1. (3.22) 
When terms up to the second derivatives are included the inequality becomes 



^(o)p'(o)*r(o) + 4fy(p)p*(p)Hf(p) + {^PFfyM 

+ (* S *PH$)?(P) + (<I> f e PFf )f (0) + faBPH*)? (0) < 1. (3.23) 



In this case the domain of the allowed values for the charge radius and the convexity is 
the interior of an ellipse in the p 2 — c plane of the form 

(p 2 - p 2 ) 2 + S[(c - c) - T(p 2 - p 2 )} 2 = K 2 , (3.24) 

where p 2 , c, S, T and K are readily determined through a comparison with eq. (3.23). 

4. RESULTS 



In this section we discuss the numerical results of our analysis. The lower and upper 
bounds on the charge radius p 2 is obtained from Eq.(3.19) with basically reliable inputs. 
Actually, most of the predicted values of the charge radius from various models are well 
below the upper bound p 2 < 4.5. The lower bound p 2 > —1.9 is not interesting, since from 



Bjor ken's sum rule ||30||, the positivity of the charge radius is expected. On the other hand, 
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Eq.(3.24) gives the correlation between the charge radius p 2 and the convexity parameter 
c. In Eq.(3.24), the numerical values of the parameters are: 

f = 1.31, S = 256, c = 1.45, T = 1.56 and K = 3.22. (4.1) 

With these parameters, the solid curve in Fig. 1 restricts the allowed range for the charge 
radius and the convexity parameter of the IW function to lie within the interior of an 
ellipse in the p 2 — c plane. 

This result can be used to test the compatibility of some phenomenological models 
for the IW function. Each model leads to definite values for the charge radius and the 
convexity. For illustration, we consider the following models. 

(1) The MIT Bag model |3l| provides the form of IW function as 

eM = (4r) 3 ' 5 ^- ( 4 - 2 ) 

LO + 1 

This model gives the values of charge radius p 2 and convexity parameter c as 2.35 and 
3.95 respectively. 



(2) In Simple Quark Model | 32fl , the IW function is well approximated by the 
formula 

= (^t) L32+ ^, (4-3) 
which predicts the computed value of the charge radius and the convexity parameter to 
be 1.0 and 1.11 respectively. 

(3) QCD Sum Rule predicts 

£H = (^)*exp[-0.8^]. (4.4) 
10 + 1 uj + 1 

The value of the charge radius p 2 = 0.65 predicted in this model, gives the corresponding 

value of the convexity to be 0.47. 

(4) The IW function in the Skyrme model |34j has the form 

f (w) = 0.99exp[-1.3(cj - 1)], (4.5) 
10 



with p 2 = 1.3 and c = 0.85. 

(5) The Relativistic Three Quark Model |35| calculation gives the approximate 



form of IW function near the zero recoil point to be 

*M = (^T) L7+ ^ (4-6) 

which indicates the numerical values p 2 = 1.35 and c = 1.75. 

(6) In Infinite Momentum Frame Quark Model fl36| , the overlap integrals for 
the hadronic wave functions of parent and daughter baryon calculated in the infinite 
momentum frame leads to the following form of the IW function 



£H = (-) 7 exp -k — — ] — v — , (4.7) 

where 

tf^) = i°° dy{y + x) l e- z \ (4.8) 

with m — 1 and «; = 1.5. This gives the corresponding values of p 2 and c to be 3.04 and 
6.81 respectively. 

In Fig. 1, the corresponding values of charge radius and convexity parameter in these 
models are plotted as points. In Table I, the predicted value of p 2 and the computed value 
of c parameter in these phenomenological models are shown. The fourth column of the 
table gives the allowed range of c for the predicted value of p 2 for each model to test the 
compatibility of the corresponding model with the dispersive bound. 

The predicted values of the charge radius and convexity parameter in the Simple Quark 
Model and QCD sum rule lie well within the allowed domain, provided by dispersive 
approach, the others lie close to the periphery or quite outside this region. 
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5. CONCLUSION 

In this paper we have related the two-point function in QCD corresponding to the 

AfyAj, pair production (evaluated at one loop) to the elastic formfactors via analyticity, 
crossing symmetry and dispersion relations. Some techniques of complex analysis , i,e. 
conformal mapping, introduction of Blaschke factor and the Schur-Caratheodory type 
inequality are employed to derive constraints on a simple parametrization of the elastic 
formfactors. In the heavy quark limit, these elastic formfactors and the formfactors for 
the semileptonic decay A& — > A C W are connected to a single IW function through heavy 
quark symmetry. As a result, we have obtained model independent dispersive bounds 
on the charge radius and the curvature of the IW function in an expansion around the 
zero recoil point u — 1. These bounds turn out to have nontrivial implications on the 
compatibility of a number of models, which we have depicted in Table I. The constraints 
can be made more stringent with the inclusion of higher derivatives , which would give 
stronger correlation between the charge radius and the convexity parameter of the IW 
function by reducing the allowed domain in p 2 — c plane. 

In previous works [^3], |24H , this dispersive approach has been applied directly to the 



formfactors for the A& — > A C W decay. But the power in that approach was somewhat 
limited due to the possible presence of be bound states below the threshold for AtA c 
production, which have not been observed experimentally and thus require inputs from 
potential models. Our approach on the other hand involves only the elastic formfactors 
and the use of experimentally known bb bound states (Ts) enables us to impose more 
stringent bounds on the IW function without having to rely much on specific models. 
This approach can also be employed to study heavy baryon lifetimes and such works are 
in progress. 
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Figure Captions 

1. The charge radius p 2 and the convexity parameter c are constrained to lie within 
the ellipse shown in figure. Values of p 2 < are excluded by Bjorken's Sum Rule. 
The points in the figure denote the numerical values of the charge radius and the 
convexity, predicted in different phenomenological models. The numbers appearing 
next to the points in the figure correspond to the models indicated by bold faced 
letters in Section 4. 
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TABLE I 



Model 


Value 
ofp 2 


Value 
of c 


Allowed range of 
c in p 2 — c plane 


MIT bag Model 


2.35 


3.95 


2.89 - 3.27 


Simple Quark Model 


1.01 


1.11 


0.79-1.19 


QCD Sum Rule 


0.65 


0.47 


0.23-0.63 


Skyrme Model 


1.30 


0.85 


1.25 - 1.64 


Relativistic Three-Quark Model 


1.35 


1.75 


1.32-1.72 


Infinite Momentum Frame Quark Model 


3.04 


6.81 


3.98-4.32 
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